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We report on polarization property discrimination of inter Landau level transitions using 
reflection- type optical Hall effect measurements from 600 to 4000 cm -1 on epitaxial graphene grown 
by thermal decomposition of silicon carbide. We observe symmetric and anti-symmetric signatures 
in our data due to unpolarized, i.e., isotropic, and polarized, i.e., anisotropic Landau level transi- 
tions, respectively. The isotropic transitions shift with square-root magnetic field dependence typical 
for stacks of decoupled graphene mono-layers. The anisotropic Landau transitions match those pre- 
viously observed in Bernal stacked bi-layer graphene, and predicted recently for Bernal stacked 
tri-layer graphene. The anisotropic signatures can be explained by coupling with free charge carrier 
magneto-optic plasma oscillations within the transition layer at the substrate graphene interface. 
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Epitaxial graphene grown by thermal decomposition 
onto SiC substrates has received tremendous interest 
due to its unique physical and electronic properties. 0- 
7J In particular, infrared magneto-optical spectroscopy 
has been widely applied to probe the electronic states of 
graphene in the vicinity of the Dirac point, the Fermi 
velocity and the mobility of free charge carriers which 
behave as quasi Dirac particles, for instance, [ll 13. lil-ficij] 
However, despite these efforts the experimental determi- 
nation of the complex- valued magneto-optic permittivity 
tensor e MO in the quantum regime which governs these 
effects has not been reported so far. The on-diagonal 
£ mo com p nents are conveniently accessible by experi- 
ment and have been successfully related to the graphene 
bandstructure. 11, [13] O n the other hand, information on 
the off-diagonal e MO components is scarce although these 
elements have been predicted to yield crucial informa- 
tion particularly when probing multilayer graphene. [l~3j| 
Faraday rotation measurements at infrared wavelengths 
currently offer the most comprehensive insight, and giant 
Faraday rotation effects and correlations with the Dirac 
Landau level structure were reported. [liHI(J Faraday ex- 
periments provide rotation and intensity, and allow for 
partial determination of £ MO only. For full access to e MO , 
required to evaluate the polarization properties of Landau 
transitions in epitaxial graphene, independent informa- 
tion for on- and off-diagonal components of e MO must be 
obtained from experiment . Sufficient information is pro- 
vided by optical Hall effect measurements, as discussed 
in this letter. 17, IH We report here on our observation of 



gins by reconstructing e MO using simple model scenarios. 

The optical Hall effect determines changes of optical 
properties of thin film samples under the influence of ex- 
ternal magnetic fields. 18l42G| In contrast to Faraday rota- 
tion, measurements are conveniently taken in reflection- 
type arrangement at oblique angle of incidence, thereby 
discriminating between parallel and perpendicular polar- 
ization. Measurements are performed in the Stokes vec- 
tor approach, and results are reported in the Mueller 
matrix presentation, which allows immediate differenti- 
ation between polarization conserving (isotropic) as well 
as polarization cross-coupling (anisotropic) sample prop- 
erties. Subsequent data analysis using model approaches 
for the dielectric function, or equivalently the optical con- 
ductivity, permits quantitative access to physical model 
parameters. In the Stokes formalism the Mueller matrix 
M connects the Stokes vector of incident and reflected 
electromagnetic waves S in and S out , respectively, where 
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isotropic and anisotropic Landau transitions from optical 
Hall effect measurements which differ in their magnetic 
field dependencies, and we discuss possible physical ori- 



The optical Hall effect determines 
magnetic field-induced differences <5M with respect to M 
at zero field. <5M contains non-zero on-diagonal block 
elements only (SM n , SM W , 5M 2 i, SM 22 , SM 33 , SM 3i , 
5M43, SM44) when the magnetic field induced sample 
response is purely isotropic (no cross-polarization), and 
additional, non-zero off-diagonal block elements (6 Mis, 
SM 31 , SMu, SMn, SM 23 , SM 32 , 8M24, SM 42 ) when the 
magnetic field induced sample response is anisotropic, 
i.e., polarization mode coupling occurs. [2l| We note that 
in our setup elements of both 4 th row and column are 
inaccessible, and all remaining elements are normalized 
by Mn removing light source base line fluctuations, pro- 
viding 8 independent information. 

The epitaxial graphene sample investigated in this let- 
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FIG. 1. Optical Hall effect ex- 
perimental data (green) and best- 
model fit data (red) for epitaxial 
graphene at B c =+(5. 66+0.02) T 
and T — 1.5 K. The angle of inci- 
dence is <P a = 45°. DTC, BLG and 
SLG denote contributions assigned 
in this letter to Drude-type carri- 
ers (DTC), bi- layer graphene LL 
(BLG) and single-layer graphene 
LL (SLG), respectively. Signatures 
indicated by SLG are isotropic and 
do not occur in off-diagonal-block 
elements SM32, SM32, SMs2, and 
8M32, while features labeled with 
DTC and BLG are anisotropic and 
cause cross-polarization. 



ter was grown by sublimation on the C-polar (0001) sur- 
face of a semiinsulating 6H-SiC substrate. During the 
growth, the SiC substrate was heated to 1400 °C in ar- 
gon atmosphere. Further information on growth condi- 
tions can be found in Ref.|22j. Optical Hall effect measure- 
ments where carried out at 4> a = 45° angle of incidence in 
the spectral range from 600 to 4000 cm -1 with spectral 
resolution of 1 cm -1 .[l7| The sample was held at tem- 
perature T = 1.5 K. The magnetic field was varied from 
B = T to 8 T in 0.1 T increments while the magnetic 
field direction was parallel to the reflected beam result- 
ing in a magnetic field B c = B/y/2 parallel to the sample 
normal. Quantitative optical Hall-effect data analysis re- 
quires stratified layer model calculations in which param- 
eterized dielectric functions are used. Least-square prin- 
ciple methods are employed in order to vary model pa- 
rameters until calculated and experimental cllipsometric 
data are matched as closely as possible (best-model). [23| 
The dielectric functions for silicon carbide are composed 
of Lorentzian-broadcned oscillators described by the re- 
spective longitudinal-optical (LO) and transverse-optical 
(TO) phonon frequencies. [24]] 

Figure [T] depicts results of optical Hall effect measure- 
ments at B c = +(5.66±0.02) T. Graphs are arranged ac- 
cording to their appearance in the Mueller matrix, with 
the top left corner (6Mu) omitted. The spectral response 
observed in the on-diagonal block elements is distinctly 
different from the off-diagonal block response. Compar- 
ing representative elements, e.g., and SM32, two 
sets of pronounced resonances which wc attribute be- 
low to Landau level transitions, can be identified. The 
first set of transitions, indicated with vertical arrows la- 
beled LL SLG , is isotropic, i.e., not associated with cross- 
polarization, and spread out over the entire measured 
spectral range. These transitions do not occur in off- 
diagonal block elements. The term "SLG" stands for 
single-layer graphene. The second set of resonances, indi- 
cated with vertical arrows labeled LL BLG , is anisotropic, 
i.e., associated with cross-polarization, and occurs in a 
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FIG. 2. Selected on- and off-diagonal-block optical Hall-effect 
spectra (T = 1.5 K, = 45°) for B c = +(0.707 ± 0.002) T 
to B c = +(5.66 ± 0.02) T in 0.07 T increments. The graphs 
are stacked by 0.006. (a) Isotropic Landau level tran- 

sitions, indicated by letters according to Sadowski et al. (b) 
5Mz2', Anisotropic Landau transitions, indicated by blue and 
red arrows. 

narrower range from 600 to 1500 cm" 1 . The term "BLG" 
indicates bilayer-layer graphene, however, as will be dis- 
cussed below, a subset of these transitions will be assigned 
to tri-layer graphene. In addition, a pronounced feature 
observed at 970 cm -1 , indicated with vertical arrows and 
labeled by DTC, is common to all graphs in Fig. [1] and 
hence anisotropic. The term "DTC" stands for Drude- 
type carriers. 

Non-zero off-diagonal block Mueller matrix elements 
are inherently tied to the existence of off-diagonal com- 
ponents in e M °. Thus, a priori, one must conclude that 
transitions LL SLG originate from processes in sample re- 
gions with polarizability contributions to e MO that are di- 
agonal and hence isotropic. Likewise, transitions LL BLG 
are to be described by contributions with non-diagonal 
components in £ M °. Without loss of generality, the op- 
tical response of electronic systems with bound and un- 
bound excitations subjected to external magnetic fields 
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can be constructed by using magneto-optic polarizability 
functions \+ an< i X- f° r right- and left-handed circularly 
polarized light, respectively, [l^ In Cartesian coordinates 
£ MO is then anti-symmetric: 

(X++X-) Kx+-x-) o\ 
-Kx+-x-) (x + + x-) o , (l) 

0/ 

where I indicates the unit matrix and the magnetic field is 
taken along the z-direction. All other dielectric contribu- 
tions are omitted. For x+ X— > £M ° describes a medium 
with anisotropic magneto-optical properties which pro- 
duce non- vanishing off-diagonal block elements in the op- 
tical Hall effect, whereas for x+ = X-i £M ° describes a 
medium with isotropic magneto-optical properties, and 
off-diagonal block elements in the optical Hall effect do 
not occur. Hence, for transitions LL SLG \+ = X-i while 
for LL BLG x+ 7^ X-- A semi-classical description for \± 
can be obtained using n series of Lorentzian-broadened 
Green functions at energies fiwo.n with spectral weight 
scattering life time l/ 7 n and including coupling to 
a magneto-optic plasma mode ui c 



X± 



2 

0,n 



uj — iujj p-n ± ILJUJc 



(2) 



When the plasma coupling is turned off, i.e. 



-> o, 

£ mo k ecomcs symmetric. Standard layer model calcula- 
tions to determine the Mueller matrix elements of epi- 
taxial graphene on silicon carbide using Eqs. ^ and (J2J 
reproduce lineshape and isotropy of all features labeled 
LL SLG in Fig. [TJ When coupling with plasma modes is 
considered, i.e., uj c > 0, features are mapped out onto the 
off-diagonal-block elements, and lineshapes match excel- 
lently with the experimental data for transitions labeled 
by LL BLG . 

Figure [2] summarizes on-diagonal and off-diagonal- 
block Mueller matrix difference spectra as a function 
of the magnetic field. The anisotropic resonance at 
970 cm -1 , labeled with DTC, increases in amplitude with 
increasing magnetic field strength. The wavenumber at 
which this resonance occurs does not vary with the exter- 
nal magnetic field strength. The physical origin of this 
resonance is the coupled motion of bound charge displace- 
ment near the longitudinal optical phonon mode of the 
silicon carbide substrate with a free charge carrier plasma 
at the sample surface producing resonant magneto-optic 
birefringence. 3, 23, 25| Eqs. Q] and [2] can be used to ren- 
der this phenomenon when hu>o is the (bound) phonon 
mode energy. The resulting polarization contribution is 
anisotropic and occurs in all Mueller matrix elements. 
We attribute this mode to highly doped graphene layers 
in close vicinity of the interface between the substrate 
and the epitaxial graphene, originating from SiC charge 
transfer. [26[ 

The magnetic field-dependent measurements reveal 
that energy spacings of transitions LL SLG scale with the 
square root of transition index n and magnetic field B, 
indicative for the Dirac-type bandstructure with Fermi 
level close to the charge neutrality point. Data in 
Fig- El (a) are parameters ujq.u obtained from best-match 



model analysis of the optical Hall effect spectra as a 
function of B c . We attribute these transition to orig- 
inate within regions of the epitaxial graphene that are 
composed of decoupled, quasi neutral graphene sheets. Q 
Energies in Fig. [3] (a) follow Ej£ a (n) = sign(n) E ^ t /\n\ 
with E = cy/2Hefj,n\B c \ and average velocity of Dirac 
fcrmions c. The naming convention used to indicate tran- 
sitions in Fig. (3Ja) is adapted from Sadowski et al.f2l\ 
Optical selection rules for transitions between levels n' 
and n require \n'\ = \n\ ± 1. The best- match model ve- 
locity obtained from matching all data in Fig. [3] (a) is 
c = (1.01 ± 0.01) x 10 6 m/ s, in very good agreement with 
Refs.iMELl l28l . andl29l The corresponding best-match 
functions E^ G versus B c are plotted as solid lines in Fig. [3] 
(a). 

Transition energy parameters obtained from best- 
match model analysis of the off-diagonal-block optical 
Hall effect data are plotted in Fig. [3] (b) and exhibit sub- 
linear behavior. The magnetic field scaling of the en- 
ergy spacings for the anisotropic transitions suggests bi- 
and tri-layer graphene as their physical origin. The tran- 
sition energies of N-layer graphene have been described 

as [13, H3, IHH 



sign(«)-| 



(A N7 ) 2 + (2|n| + l)£ 2 



1/2 



+^ (A n7 ) 4 + 2 (2 M + 1) £ 2 (A n7 ) 2 + 



F 4 



(3) 



with coupling constant 7, layer parameter A N [31| and 
where /i = — 1 , + 1 corresponds to the higher and 
lower subbands in the limit of zero magnetic field, 
respectively. [3l[ Optical selection rules are the same as 
for E£ G (n). Using c = (1.01 ± 0.01) x 10 6 m/s, best- 
match model functions are plotted in Fig. EJb) for N = 2 
and N = 3 as blue solid and red dashed lines, respec- 
tively. No transition was observed that would correspond 
to N > 3. These transitions can be assigned to Bernal 
stacked bi-layer graphene (N = 2) and tri-layer graphene 
(N = 3). The best-match model parameters obtained 
here are 7 = (3120 ± 175) cm -1 for N = 2, corrobo- 
rating the result obtained by Orlita et al. from FTIR 
transmission experiments . j 1 Oj For N = 3, we observe 
7 = (3150 ± 20) cm -1 , which render the first experimen- 
tal confirmation of the theoretical predictions by Koshino 
and Ando for tri-layer graphene. [31( We note that these 
anisotropic Landau transitions are only observed and re- 
solved for n = 2 . . . 6 for bi-layer graphene, and for n = 4 
and n = 5 for tri-layer graphene. At this point, we do 
not know why transitions with n > 6 for N — 2 and 
n 7^ 4, n 7^ 5 for N = 3 cannot be observed. However, 
the fact that these transitions appear with anisotropic 
optical Hall effect signatures suggests their coupling with 
free charge carriers within the sample. We propose that 
stacking order defects within the bi-layer and tri-layer 
graphene allow for coupling of Dirac particles within their 
Landau levels with the cyclotron resonance of the free 
carrier plasma resulting in anisotropic Landau transition 
signatures. 

In conclusion, we determine for the first time the full 
polarization properties of Landau transitions in epitax- 
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FIG. 3. Symmetric (a; isotropic) and anti-symmetric (b; 
anisotropic) Landau transition energies in epitaxial graphene 
determined from optical Hall effect measurements at 1.5 K. 
Solid lines denote single-layer graphene (a; black), bi-layer 
(b; blue) and tri-layer graphene (b; red) dependencies on the 
magnetic field. 

ial graphene using the optical Hall effect. Magnetic- field 
dependent measurements revealed that isotropic tran- 
sitions originate from decoupled graphene mono-layers, 
while anisotropic transitions reflect dependencies previ- 
ously recorded for bi-layer graphene. Furthermore, we 
observe anisotropic transitions previously predicted for 
tri-layer graphene. The anisotropic transitions corre- 
spond to optical conductivity contributions with non- 
vanishing off-diagonal tensor components. We suggest 
that Landau transitions, while intrinsically isotropic, pro- 
duce anisotropic conductivity contributions when cou- 
pling occurs between Landau levels and free charge car- 
rier modes. We suggest that such coupling may be 
promoted via stacking defects in bi-layer and tri-layer 
graphene. The knowledge of the polarization dependent 
optical response of Landau-level transitions is found in- 
valuable for accurate assignment of the origin of the Lan- 
dau transitions occurring in epitaxial graphene. 
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